Fluctuations of the luminosity distance 
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We derive an expression for the luminosity distance in a perturbed Friedmann universe. We define 
the correlation function and the power spectrum of the luminosity distance fluctuations and express 
them in terms of the initial spectrum of the Bardeen potential. We present semi-analytical results 
for the case of a pure CDM (cold dark matter) universe. We argue that the luminosity distance 
power spectrum represents a new observational tool which can be used to determine cosmological 
parameters. In addition, our results shed some light into the debate whether second order small 
scale fluctuations can mimic an accelerating universe. 
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I. INTRODUCTION 

Some years ago, to the biggest surprise for the physics 
community, measurements of luminosity distances to far 
away type la supernovae have indicated that the Uni- 
verse presently undergoes a phase of accelerated expan- 
sion If the Universe is homogeneous and isotropic, 
i.e. a Friedmann-Lemaitre universe, this means that the 
energy density is dominated by some exotic 'dark energy' 
which obeys an equation of state of the form P < — p/3. 
The best known dark energy candidate is vacuum energy 
or, equivalently, a cosmological constant. This discov- 
ery has lately been supported by several other combined 
data sets, like the cosmic microwave background (CMB) 
anisotropics combined with either large scale structure or 
measurements of the Hubble parameter 0]. 

On the other hand, since quite some time, it is known 
that locally measured cosmological parameters like Hq or 
the deceleration parameter qo might not be the ones of 
the underlying Friedmann universe, but they might be 
dressed by local fluctuations |3(. Therefore, it is of great 
importance to derive a general formula of the luminosity 
distance in a universe with perturbations. To some ex- 
tent, this has been done in several papers before 0,13. 
But the formula which we derive here is new. We shall 
comment on the relations later on. 

Lately, it has even been argued that second order per- 
turbations might be responsible for the observed acceler- 
ation and that no cosmological constant or dark energy 
is needed 0,0- This claim is very surprising, as it seems 
to require that back reaction leads to big perturbations 
out to very large scales, contrary to what is observed in 
the CMB. This proposal has thus promptly initiated a 
heated debate 

On the one hand, the present work is a contribution 
in this context. We calculate the measurable luminosity 



distance in a perturbed Friedmann universe and deter- 
mine its fluctuations (within linear perturbation theory) . 
We show that these remain smaller than one and there- 
fore higher order perturbations are probably not relevant. 
The main point of our procedure is that we use only mea- 
surable quantities and not some abstract averaged expan- 
sion rate to determine the deceleration parameter. We 
actually calculate the luminosity distance dz,(n, z) where 
n defines the direction of the observed supernova and 
z its redshift. We then determine the power spectrum 
Ci(z, z') defined by 



d L (n,z) = y j a im {z)Yt m {n) 



( III 



Ci(z,z') = (aim(z)a} m (z')) 



(1) 
(2) 
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Here the (•) denotes a statistical average. Like for the cos- 
mic microwave background, statistical isotropy implies 
that the C/s are independent of m. 

We then analyze whether the deviations of the angu- 
lar diameter distance from its background value can be 
sufficient to fake an accelerating universe. 

Aside from this problem, the new variable which is 
defined and calculated in this paper, might in principle 
present an interesting and novel observational tool to de- 
termine cosmological parameters. And this is actually 
the main point of our work. We hope to initiate a new 
observational effort, the measurement of the luminosity 
distance power spectrum, with this paper. A detailed 
numerical calculation of the cIl power spectrum and the 
implementation of a parameter search algorithm are post- 
poned to future work. Here we simply show that for large 
redshifts, z > 0.4 and sufficiently high multipoles, I > 10 
the lensing effect dominates. However, at smaller redshift 
and especially at low ts other terms can become impor- 
tant, most notably the Doppler term due to the peculiar 
motion of the supernova. 

The paper is organized as follows. In Section [H] we 
derive a general formula for the luminosity distance valid 
in (nearly) arbitrary geometries. In the next section we 
apply the formula to a perturbed Friedmann universe. In 
Section lTVI we derive general expressions for the g?l power 
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spectrum in terms of the Bardeen potentials. We then 
evaluate our expressions in terms of relatively crude ap- 
proximations and some numerical calculations for a sim- 
ple tt M = 1 CDM model in Section El In Section IVT1 wc 
discuss our results and conclude. 

Notation: We denote 4-vectors by arbitrary letters, 
sometimes with and sometimes without Greek indices, 
k = (k 1 - 1 ). Three-dimensional vectors are denoted bold 
face or with Latin indices, y = (y l ). We use the metric 
signature (—,+,+,+). The covariant derivative of the 
4-vector k in direction of the 4-vector n is often denoted 
by V n k ee Kfc%). 

II. THE LUMINOSITY DISTANCE IN 
INHOMOGENEOUS GEOMETRIES 




FIG. 1: A light beam emitted at the source event S ending 
on the observer O. At the source position, the plane normal 
to the source four-velocity is indicated. 

We consider an inhomogeneous and anisotropic uni- 
verse with geometry ds 2 — g^dx^dx" . We place a stan- 
dard candle emitting with total luminosity L ( energy 
per unit proper time) at spacetime position S. Its four- 
velocity is us- An observer at spacetime position O with 
four velocity uo (see Fig. ^| receives the energy flux F 
(energy per unit proper time and per surface) . The lumi- 
nosity distance between the source at S and the observer 
at O is defined by 



d L (S,0) = 




The observer measures the flux F and 'knows' the intrin- 
sic luminosity L of the standard candle. Furthermore, 
she determines the source redshift z and direction n and 
thereby obtains the function d^ (n, z) , which we now want 
to express in terms of the spacetime geometry. 

Be dQs the infinitesimal solid angle around the source 
and dA{x) the infinitesimal surface element on the sur- 
face normal to the photon beam at the position x along 
the photon trajectory from S to O, then 

d K S > °) = + z ? = l det J (°> S )K 1 + *? ■ ( 4 ) 



Here J is the so called Jacobi map mapping initial direc- 
tions 59g around the source into vectors 5x Q transversal 
to the photon beam at the observer position [jj , 

8x> = J» 0l {0,S)86% . (5) 

The factor 1 + z = ujs /ojo is the redshift of the source. 
There is a factor 1 + z due to the redshift of the emitted 
energy and a second factor due to the time dilatation in 
F oc dEo/dro with respect to L = dEs/drs- If k denotes 
the 4-vector of the photon momentum and us and u<j are 
the source and observer 4- velocities respectively, we have 

- us = (k ■ u s ) = gflv (S)k»(S)u s (S) and (6) 

-a Jo = (k-u o )=g^(0)k»(0)u o (0) (7) 

If we have a standard candle source of which we 
know L and we measure F, we can therefore deter- 

1/2 

mine |det J(0, S)\ ujs/ojo , which contains information 
about the spacetime geometry. Of course it also depends 
on the source and observer velocities. The Jacobi map 
J^ a (0, S) maps direction vectors normal to the photons 
direction and normal to us at S into vectors normal to 
the photon direction and uo at O. It depends on the 
source velocity us and on the curvature tensor along the 
photon geodesic from S to O. As we shall see, it does 
not depend on the observer velocity uo- 

Even though in the form JSJ, J is given by the 4x4 
matrix J7 M Q (0, S), we have to take into account that the 
vectors Sxq as well as 89g live in the two dimensional 
subspace normal to uo respectively us and normal to 
the photon direction at O and S. The latter are given by 

no = — (UO) + (k(0) ■ u )u ) and (8) 

UJ 

n s = — (k(S) + (k(S) ■ u s )us) . (9) 

UJ S 

The photon direction vectors ns and no are normalized 
spacelike vectors pointing into the photon direction in the 
reference frame of the source at S and of the observer at 
O respectively. Denoting the projectors onto the sub- 
spaces normal to us,ns and uo,no by Ps and Po we 
have 

{PsYv = ^v+u^usv-n^nsv and (10) 
{PoYu = ^ v + u u 0l/ -ri n 0v . (11) 

The true Jacobi map is J(0, S) = PqJPs understood 
as two dimensional linear map. For convenience we shall 
write it as four-dimensional application and determine 
its determinant as the product of the two non-vanishing 
eigen- values. 

To determine the Jacobi map we now derive a differ- 
ential equation for the evolution of the difference vector 
Sx fl (X) in a given direction 56 $ along the photon trajec- 
tory. The final value Sx^^Xo) then depends linearly on 
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the initial conditions 



For this we denote the pho- 



ton trajectory by /"(A, 0) and parameterize neighboring 
light-like geodesies by f a (X,5y). The 4- vector 



k a (Sy) 



df a (\,6y) 
dX 



is the tangent of neighboring photons at Sy and 

df a 
Sx a = -^—5y l 

connects the geodesies / a (A, 0) and / a (A,<5y). Since the 
'beam' f a (X,y) describes photons which are all emitted 
at the same event S they have the same phase (eikonal) 
S. With k a = — V Q <S we therefore have 



= V Sx S = Sx a V a S = -Sx a k a 



(12) 



In order for the 4- vectors Sx a (y) to sweep a surface nor- 
mal to uo at the observer event O at A = Xo, we also 
need (dx(Xo) ■ uo) = 0. This is a priori not true. How- 
ever, we can re-parameterize / by 



A -> A = A + h(y) and y -> y = g(y) 



(13) 



Under this reparameterization Sx transforms as 5x a — > 
Sx — 5x a + k a Sh. It is easy to see that g a p5x a 8x" = 

gapbx Sx^ \ hence the length of the vector Sx is invari- 
ant under this reparameterization. Since uo is timclikc, 
(fc(Ao) • uo) ^ and we can hence choose a parameteri- 
zation such that (Sx(Xo) ■ uo) = 0. 
The directions 59° are given by 

SO" = — (V k 5x) a = —{V Sx k) a . (14) 
UJ S uj s 

The last equality requires a brief calculation which can 
be found, e.g. in Q. To convince oneself that the above 
definition of S9 a is suitable, one easily verifies (see Q) 



that 



is normal to the source velocity us and the 



photon direction ns and that it is normalized. 

To find the differential equation for 5x(X) we use the 
relations 



i?V„fc' 3 - (V M V„ - V„V„) k a 
(W k Sx) f3 = fc Q V Q< 5x' 3 = 5x a W a k f3 = (W Sx kf .(15) 

Furthermore, 

R% v k k»5x u = W8x v (V„V„ - V„V„) k a 

= Sx u V k (V,fc a ) - k»V 5x (V M fc Q ) 
= V fc (5x»V„k a ) - {V k 5x v ) {V u k a ) 

-V 5x (fc^fc") + (V Sx k») (V M fe a ) 
= S7 k {8x v S7 v h a ) 

= Vfc {k u V v Sx a ) = V k (u> s S9 a ) . (16) 



From the third to the fourth line we have used that 
Vfcfa = Vg x k and V k k = 0. We therefore obtain the 
system of equations 

V k (uj s 59 a ) = Rfj^WSx" (17) 
V k (Sx a ) = oj s 59 a . (18) 

With the definition of the covariant derivative this fi- 
nally gives 



^1 = -T« v mx» + u> s 86° 
= C*{X)8x u +uj s 69 a 

d(uj s S9 a ) 



(19) 



dX 



RS^tfWx" - T°k"Lj s 60 v 



= A^{X)5x v + C^{X)u s 86 v (20) 
where we have set 

C£(A) = -I^fc" and A a p {\) = B% tlft WW ■ (21) 



We now define 



Z = 



( Sx c 



\ uj s 59' 

This (8 component) vector then satisfies the equation 



(22) 



dZ(X) 
dX 



= B(X)Z(X) 



with 



~ ' A%X) C£(A) 



(23) 



(24) 



The initial conditions are Sx a (Xs) = since all photons 
start from the same source event and (k a 69 a ) (As) = 
(tig S9 a (Xs j) — as we have seen above. The solution of 
Eq. ill'MI) therefore provides a linear relation between the 
initial condition S9 a (Xs) and 6x a (X), 



6x a (X) = J?(X)56P(X s ) 



(25) 



With J(Xo) we can then easily determine the true Jacobi 
map J(0,S) = P J{X )Ps- 



III. THE LUMINOSITY DISTANCE IN A 
PERTURBED FRIEDMANN UNIVERSE 

A. Conformally related luminosity distances 

We consider two geometries related by 

ds 2 = g^ v dx^dx v — a 2 (x)g flv dx fJ 'dx l/ = a 2 (x)ds 2 . (26) 

We want to relate the angular diameter distances of the 
two metrics. If k is a light-like geodesic for the metric 
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ds with affine parameter A, then k — ark is a light-like 
geodesic for ds 2 with affine parameter A determined by 

dX 2 

Tx =a ■ 



Furthermore, be u p 



df 



the 4-velocity of an observer 



with metric ds 2 and be f its proper time such that 
g l _ LU u^u v = —1, then = #r is the corresponding 4- 
vector of the observer with respect to the metric ds 2 with 
proper time r if 4^ = a. In other words 



dx^ dx* dr 



df 



dr df 



(27) 



The redshift of a photon emitted at S and observed at O 
determined with respect to the two metrics is therefore 
related by 



1+5 



wo (g^k^u u )o asig^ukf-u^o a s 



(1+z) . 
(28) 

To determine the relation between the Jacobi maps 



J" 



we just have to remember that angles are 

not affected by conformal transformations, but distances 
scale with the conformal factor a. Therefore 



fif a 



a °T^3 = a oJ(S,0) 
S9 s 



(29) 



det J(S,0) = a 2 o detJ(S,0) . (30) 
For the angular distance relation we finally obtain 



d L = (l + f)y|det J(S,0)\ 

2 2 

= ^(l + z) v /|det J(S,0)\ = ^d L . (31) 
as as 

This relation is very useful in Friedmann cosmology. 
The Friedmann metric is given by 



ds 
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(—di] 2 + "fijdx l dx^ = a 2 ds 2 



(32) 



where 7 is the metric of a 3-space with constant curva- 
ture K. The luminosity distance of a photon emitted at 
conformal time 775 and observed at r\o with respect to the 
metric ds 2 is simply r\o ~ Vs = f V ° drj. The Friedmann 
equation for a universe containing matter, radiation, cur- 
vature and a cosmological constant reads 



— I — Hn \Q m a 1 



n iad a- 2 + n K + n A a 2 ] , (33) 



where we have normalized ao = 1 and we have intro- 
duced the density parameters fi m = p m (vo) I Pc(vo) , 
^rad = Pr a d(r]o)/Pc(r]o), = - K '/ Hj 2 and SI a = 

A/(3H 2 ). 

After the variable transformation to z + 1 = 1/a, dz — 
—da I a 2 we obtain 



drj = 



H^ 1 dz 



^a-ad(i + zY + n m {z + 1) 3 + n K (z + 1) 2 + n A 



This leads to the well known expression for the luminosity 
distance to an object emitting at redshift zs observed 
today at zq = 0, 



dbizs) 
1 + zs 



Friedman _ ^0 VS 



as 



Z S +1 



dx 



Hq J l V^rad^ 4 + fl m X 3 + Qk% 2 + ^/ 



(34) 



Comparing this expression with the measured luminosity 
distance from supernovae type la at different redshifts 
has led to the claim that the cosmological constant be 
non- vanishing 0. 



B. The Jacobi map in a perturbed Friedmann 
universe 

We now consider a Friedmann universe with scalar per- 
turbations. In longitudinal (or Newtonian) gauge the 
metric is given by 

g^dx^dx" = a 2 [-(1 + 2^)dr] 2 + (1 - 2$)~fi j dx i dx j ] . 

(35) 

For perfect fluids the metric perturbations '5 and $ are 
equal. We assume in the sequel $ = \P. Furthermore, 
we consider a spatially flat universe (K = 0), so that 
Jij=Sij. 

We now determine the luminosity distance for the met- 
ric 

ds 2 = -(1 + 2#)<V + (1 - 2*)<%cfa:W . (36) 

We then relate this to the physical luminosity distance 
via the relation 13111 . 

We assume that the galaxy containing the supernova as 
well as the one containing the observer are moving with 
the cosmic fluid. To first order in the perturbations, the 
four velocity of the cosmic fluid is given by 

(u") = 0—*,v i ), (37) 

where v l is the peculiar velocity field. 

1 . Redshift 

The photon geodesic is obtained by integrating the 
geodesic equation to first order. Since the background 
is Minkowski, the background photon momentum is con- 
stant and we may normalize the affine parameter such 
that k° = 1 and k l — n l with n 4 n 4 = 1. Here over- 
bars denote background quantities. For the perturbed 
4-velocity of the photon we may still assume kg = 1. 
The geodesic equation then gives (to first order) 

fc (A o ) - k°(X s ) = fc (A o ) - 1 = -2 f ° dAVtf (A) • n 
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rXo 

= - 2*|° + 2 / d\ii and 

fc l (A G ) - fc*(A s ) = 2n l (* - * s ) - 2 / ° dAfl^A) . 

The redshift of a photon emitted at spacetime position 
S and observed at O then becomes 

(38) 



4* — —4° 

(40) 

The Christoffcl symbols and the Ricci tensor of the 
perturbed metric are given in Appendix El Spatial in- 
dices i or j are raised and lowered with the flat metric Sij . 
Therefore, no special attention is paid to their position. 

To solve it, we now split the system (12. '{I) into its zeroth 
and first order components, 



2. The perturbed Jacobi map 



Z = Z (0) + Z {1) and B = B + B {1) (41) 



To determine the Jacobi map we have to solve the sys- 
tem (I23[l to first order. We first determine the maps C a p 



and A a p which make up the matrix B 



N 



M ■ 



According 



to Eq. C a p = -r^jfeT. Since is already first 

order, we may insert the zeroth order expression for k 1 
leading to 



C? = 

a = 



(39) 



Here we denote the derivative along the geodesic with 
a prime and the derivative w.r.t. conformal time by an 
over-dot, ^ = ' and J^- = '. The matrix A is given 
by Ap 1 = R^gk^k". Again, since R°^ v p is of first order, 
we may insert the zeroth order expression for the pho- 
ton velocity. Note that A a p, unlike C a p, is symmetric. 
Computing the Riemann tensor of our perturbed metric 
we obtain 



4° 



A» = 2 di *-^- d ± 



dX 



dX 



To zeroth order, the photons move along straight lines 
and the energy is not redshifted so that we simply obtain 
S9 a (X) = 69%, w(A) = lo s and Sx a (X) = (A - A s )u; S (50g. 
For the Jacobi map this implies Jp = (Ao — Xs)cosSp. 
The projector onto the tangent space normal to the ob- 
server velocity and the photon direction is simply P$ = 
Po = P, where 



p0 _ 

M) — 

P l = 
j 



(42) 



The zeroth order 2-dimensional Jacobi map is therefore 
given by J| = {PJP) a 



t° — t° — P — n 

j — ,J i — j — u 

Jj = (Ao - A s Vs (5) - n l n 3 ) 



(43) 



The 2-dimensional determinant of the Jacobi map is 
therefore det J = (Ao — As) 2 w|, leading to the flat space 
luminosity distance dr, — Ao — As = rjo — rjs- For the 
last equality we have used that n° = -5? = 1. In an 
unperturbed Friedmann universe this reproduces (|34l) . 

Since C and A are already first order, the first order 
differential equation becomes 



^5x a{ - 1 \X) = C" {1 \x)5x fj {X) + {oj s 59 a )^{X) 
dX 



±(oj s Se a )^(X) = A a g {1 \x)8xf ) {X)+C a B (1 \x)uj s 8¥{X) . 



dX 



(44) 



Making use of the background solution we obtain 

cX 



(S9 a ) {1) {X) 
Sx a ^{X) 



dX' (A|(A')(A' - As) + C£(A')) 6% + (S6 a s )W 



A r X f\' 

Ifioit \ l\t\l 



dA'Cf (A')(A' - As) + / dX' / dX" {A1(X")(X" - X s ) + C$(X")) 

As •' A.s •' As 

+(A-As)(^s^) (1) • 



Ld s 89 f i 



(45) 



(46) 
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The first order contribution to the unprojected Jacobi map then becomes 



cjfj^^iXo) = d\C%(\)(\-\ s )+ [°dX / A dA'(^(A')(A'-A s ) + ^(A')) . (47) 



Xs J \s J Xs 



We want to calculate 



J (1) = (PoJPs) {1) = PoJ {1) Ps + P ( o ] JPs + PoJP ( s ] • (48) 
A short calculation, inserting our results for C and A gives 

(PoJ (1) Ps) % j = U • (5) - 71%) + Wj - nWWkj - n 3 n k W l k + n l n 3 n k n l Wu (49) 

with 

U = -2* S (A - A s ) + 2 j ° dA*(A) and W tj = -2 f ° dX f dA'S^tf (A')(A' - A s ) . (50) 

J Xs J Xs J Xs 

Implicit summation over repeated (spatial) indices is assumed and n' 1 — rij, Wj = Wij = W % K 
Calculating also the first order contributions to the projections we finally obtain 

A = o 

J° = cj s (X - As) {vq - n l n k v^) 
J l = uj s (X -X s )(-v l s + n t n k v k s ) 



J) = u; s (X - A s )| ^1 - 2v&s + ^ 2 ^ £° dA*(A)^J 5) + n'rij ^ - 1 + 2* s - — 2 — dXV(X) 

-n(v + v s ) - 2 f ° dAV*(A)n + 2n • j + n l v j + n 3 v% + 2 f ° dAfit,-* (A)n* - n^g - n^^' 

I " d\ [ d\'(\' - As) (didjV - nWdjdk* - v?n k d l d k ^ + n i n j n k n l d k d l V} (A') 1 . (51) 
Jx s JXs \ J J 



Ao - As 

Like in the unperturbed case, the two eigenvalues of the Jacobi map are equal. This is due to the fact that the 
shear contribution to the Jacobi map still vanishes in first order. A short computation gives the eigenvalues a, 

r-Ao 



a - u s (X - As) j 1 - 2* s + Xq 2 ^ d\^(X) 

~ A -A S fj dX !x dX ' {X ' ~ A s)( v2vJ, ( A ') - did^^XnP) j . (52) 

The luminosity distance of the perturbed Minkowski spacetime it given by d^ — (u>s/u>o)oi- Inserting the above 
expressions and taking into account the perturbation of the emission frequency, los = ~{g^,vk iJ 'u u )s — Cos \ we 
obtain 

{2 fVo rVo 

1 - o + n ■ (v - 2v s ) + / drff + 2 / drjVV ■ n 

Vo - Vs J ns J ns 

2 f va f ri 1 f' v " f' v / \ 1 

+ / dri drj'VV-n / dr] drj'M - r? s )(V 2 * - nVM't } . (53) 

VO - VS Jns Jns VO ~ VS Jr,s Jns V 1 \ 
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Here we have also transformed the parameter A into the conformal time 77 via the relation 



dr) 
dX 



ra° (A) = 1 - 2 [ dA'V* • n 

JXs 



Now 77 is parametrizing the unperturbed photon geodesic and we interpret the potential as a function of 77, *5f(rj) — 

^(77, x(r/)). We use the notation VP = d v ^, so that = 5* + n • V*. We now also take into account expansion, which 

~ 2 
gives dh = ^-dh- 

& ^ as n 

Furthermore, we relate the peculiar velocities to the Bardeen potential via the first oder perturbations of Einstein's 
equations. Setting (u^) = a _1 (l — *,u l ) gives |Toj . 



Art) 



+ ft* 



47rGa 2 (p + p) V« 

With this we find the following result for the luminosity distance in an perturbed Friedmann universe 



(54) 



di (575,11) = — (770 - i]s){ 1 - *o + v • n 



as 



?70 - VS Jris 



no 



dr/* + 2n ■ 



1 



?7o - »?s 



no 



rV 

dr] / dTy'V* 



1 



47rGa|(p+7j)(r7 S ) V 



(WV# + V*) (77s) 



1 f'° P / 

/ dr) dr({r( -77 S )fv 2 * -iiV^jf 

'/S J ns ri.s 



(55) 



where we have introduced 7i = d/a = a = -Ha. In what follows, we further simplify the formulas by normalizing 
the scale factor to ao = 1 • 

Here we have used the linear perturbation theory solution for the source velocity vs- One might argue that the 
supernovac are highly non-linear objects inside galaxies and do not move with the velocity obtained from linear 
perturbation theory. However, we shall be interested in distances and angles which are sufficiently large so that the 
non-linear contributions to the supernova velocities are uncorrelated and therefore considering only the linear part of 
it in the correlation function is sufficient. 

Eq. ll55[l is the luminosity distance of a source in direction — n at conformal time rjs- However, this quantity is not 
directly measurable. What we do measure instead is the redshift of the source, z$ = z$ + 5z$, where z<j + l = l/a(rjs). 
Now 



d L (vs,n) = d L (r](zs),n) = d L (zs,n) = d L (z s ,n) - -^-d L (zs,n)6z s . 

dz s 



Furthermore, 



-7=-J L (z s ,n) 

dz s 



(1 + z s ) 1 d~L + Ug 1 + first order 



8z s = (1 + z s )Sz s = (1 + z s ) 
Inserting this in Eq. I|55|) leads to 



and 

rvo 

1 17s 



rvo 

*5 ~ *o + 2 / df]n ■ V* + (v - v s ) ■ n 

J vs 



(56) 



(57) 



d L (z s ,n) 



{1 f vc 
(Vo ~^)+^ (* o - v • n) - (770 - 775 + H s l ) * s + 2 J 



drj^ 



2n- 



1 rvo 



Vs 

VO rV 

dr/ / dr)'(r)' 
vs Jvs 



vo rv 



Vs J Vs 



■no - -Us 1 



Vs 



)(v 2 * - didjVr 



(58) 
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After several integrations by part, one can also derive the following expression for the luminosity distance, which can 
also be found elsewhere 0,0, where it has been derived using the evolution equations of the expansion and the shear. 



d L (zs 1 n) = {l+zs)(r)o-r)s)< [ I-7 v£T v °" n_ 1_ 7 vTT l vs ' n 

{ (vo~Vs)K s V [Vo ~ Vsjris / 

■2 I * s + ( 1 - — — ) *o 



(vo-Vs)'HsJ V (vo-vs)'H 



is 



(»7o - »7s) ^ s (»7o ~ ?7s)fts ,/„ s fao - Vs) J vs (vo - Vs) 

' d7? (^MM^)y^\ , (59) 



O70 - »7s) 



A detailed derivation of this result starting from 
Eq. H58fl is given in Appendix \B\ In this equation the 
first line, apart from the background contribution, con- 
tains the terms due to peculiar motion of the observer 
and emitter (Doppler terms). The second line can be 
identified as 'gravitational redshift'. This is, however, 
not entirely correct since this term does not vanish even 
if ^5 = \&o- The third line collects integrated effects 
proportional to line of sight integrals of "J and its time 
derivative, and the fourth and last line represents the 
lensing term with V 2 ^ cx 5 p. This term has been dis- 
cussed in the literature before • An equivalent of the 
above formula can also be found in [T^ . 

Eqs. (|58|l and (|59|l are the final expressions for the 
luminosity distance in a perturbed Friedmann universe, 
as a function of the measured source redshift z$ and its 
direction — n. In the next section we determine the lu- 
minosity distance power spectrum which is, in principle, 
an observable quantity. 
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FIG. 2: We show the dipole amplitude in a pure CDM uni- 
verse in units of the CMB dipole as a function of z for z' = 0.1, 
0.5, 1, 2 and 4 from top to bottom. 



IV. THE LUMINOSITY DISTANCE POWER 
SPECTRUM 

We now want to determine the power spectrum of the 
perturbed luminosity distance, as defined in the intro- 
duction. For notational simplicity, we drop the ~ and 
use d,L to denote the luminosity distance in a perturbed 
Friedman universe. From Eqs. (JTJ and J5J and the ad- 
dition theorem for spherical harmonics, one obtains the 
correlation function 

dh Os ) " 1 dh {zs' ) ~ 1 {dh (zs , n) d L (z s > n') ) = 
x 21+1 



C e (z 3 ,zs')Pe(n-n') (60) 



where Pi is the Legendre polynomial of order i. 



A. The dipole 

Let us first briefly look at the dipole coming from the 
peculiar motion of the observer, the term containing the 



scalar product n • vo . 
is given by 



The power spectrum of this term 



{d^(z s ,n)d^(z S ',n')) 

(z s + l)(z s , + 1) / 
\ v 



>>(n-n') 



(61) 



We assume that, like for the anisotropies in the cosmic 
microwave background, this term completely dominates 
the dipole. The luminosity distance dipole therefore has 
the same direction as the CMB dipole. To determine its 
amplitude we insert dL{r]s) = (zs + l)(vo — Vs)- We then 
obtain 



4-7T 



nj — 1 nj — 1 



(vo - Vs)(vo - Vs>) 



(62) 



The CMB dipole is given by the expression in square 
brackets. In a pure CDM universe with H = 2/rj and 
Vo/vs — V z s + 1 we obtain for the amplitude of the 
luminosity distance dipole 



Ci(z,z') =C\ 



CMB 



1 



4(V7+T- 1)(VF+T- 1) 



(63) 
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FIG. 3: We show the dipole amplitude in units of the CMB 
dipole as a function of z — z in a pure CDM universe. 



In Fig. |2 the relative amplitude of C\ as a function of 
z for different values of z' is shown. In Fig. |3| we plot 
Cx (z,z). It seems to be most promising to measure the 
dipole at relatively low redshift. But of course, the red- 
shift must be sufficiently high such that the peculiar ve- 
locities of the supernovae themselves are not strongly 
correlated with each other or with our peculiar motion. 
Hence the distance of a supernova at z or z' should be 
sufficient for linear perturbation theory to apply. This 

is safely achieved for z, z' ~ 0.1. At z = z' = 0.1 we 
have C\ (0.1, 0.1) ~ 105 x Cf MB , hence an enhancement 



of about a factor 100 with respect to the CMB dipole. 
This factor is even somewhat larger, in a A-dominated 
cosmology. Through its dependence on TL{z), measuring 
the amplitude of this dipole alone can already lead to 
new observational constraints on the expansion history 
of the universe. 



B. The higher multipoles 

We now want to express the higher C/s in terms of the 
power spectrum for the Bardeen potential. We define the 
Fourier transform 



dxe lkx <F(77,x) 



(64) 



We split the deterministic time evolution into a 'transfer 
function' T k {rj), such that #(77, k) = T k (rj)V(k). We nor- 
malize the transfer function such that linifc_>o T k (j]o) = 1. 
The power spectrum P-% of W(k) is defined by 



A: 3 (*(k)tf*(k')) = (27r) 3 £ 3 (k-k')P*(fc) 



(65) 



The <5 3 - function is a consequence of statistical homo- 
geneity. We need to determine the correlation func- 
tion of 'F for the positions x = xq — n(r)o — if) and 
x' = xo — n'(ryo — i)'). In terms of the power spectrum 
the correlation function of \& and of its derivatives as 
they enter in Eq. I|58|l can be written as (for details see 
Appendices [U] and [DJ) • 



(*(r?,x)*(r/,x') 

(n ■ V*(?7,x)*(?7 / ,x / ) 

C^\z lZ ' 
(n i n i a l d J *(ri,x)W(rf,xr) 



d nndd *\z,z' 



(V 2 *(r/,x)v]/(7/,x') 

Cf m \z,z' 
(n • W (77, x)n' • V*(r/, x') 



Ci nd * nd) (z,z' 



= y^±lcf(z, Z y,(n-n') with 

9 r rib 

= - J -jn^TMP^miKrio ~ V))je(k(r, - rf)) 
= Y?l±lrf nd *\z,z')Pe(n-n') with 

I 

= -- [ dkT k (r,)T k (r,')P^(k)j' e (k(vo-v))MKvo-v')) 



= ^±lc^ nndd9 \z, z')P e (n ■ n') with 



dkkT k ( V )T k {Ti')P 9 (k)j'ttk(Tio ~ v))MHvo - rf)) 



= E^^^Wn-n') with 



E 



dkkT k (r})T k (ri')Py(k)je(k(i)o - v))M k (Vo - v')) 



H + 1 ^nd^nd) , 1 



4tt 



(z,z')Pe(n- n') with 



= - / dkkT k (ji)T h (Tf)P 9 (k)j't(k(Tio ~ vMiHvo - rf)) 

IT 



(66) 



(67) 



(68) 



(69) 



(70) 
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(n ■ V*(77,x)V 2 *(r/',x' 



(rid*dd) 



(z,z 



(n • V* (r),-x.)n l n J didj^ \rf , x' 



^{ndty ndnd) t 



(n i rv i did j $(T}, K)n' i n' j d l d j ^(r ] ', x' 



(n7idd^> nndd) 



(z,z 



iV<9 l <9 i <i'(7 ?) x)V 2 1'(7/,x' 



C, 



(ndnd^dd) 



(z,z 



(V 2 *(t7,x)V 2 #(?/,x' 



E 



E 



2^+1 

47T 



c 



(nd>S>dd) 



(z,z')Pt(n-ri) with 



&fc 2 r fe (r?)T fc (7 ? ')P*(fc)jKfc('7o ~ v))MKvo - rf)) 



2£+l 
An 



-C 



(nd^ndnd) 



(z,z')P e (n-ri) with 



E 



2^ + 1 n (nndd^nndd) , i 



47T 



(z,z')P^(n-n') with 



- / &fc 3 r fe (^)T fc (r ? ')P*(fc)j7(fc(r/ - vM'iHvo ~ rf)) 

IT 



E 



E 



21+1 

47T 



(vidnd'fydd) / / 



(2, z')A(n • n') with 



dfcfc 3 r fc (r ? )T fe (r,')P*(fc)j7(fc(r /0 - r]))h{k(vo - rf)) 



^Cf^fvV^n.n') with 



dkk 3 T k (rj)T k (r)')P 9 (k)je(k(r)o - v))je{k(r]o - rf)) 



Using these definitions we can write the correlation function of the luminosity distance as 



(71) 



= — / dkk 2 T k {rj)T k {r 1 , )Pv{k)j[{k{r 1 o - vM(Hvo - tf)) (72) 



(73) 



(74) 



(d L (zs,n)d L (z s >,ri)} 
d L (z s )d L (z s >) 



E 



2^+1 

47T 



P e (nn>)(cf)+Cf+Cf+Cf 



a 



(75) 



where C\ collects all the contributions to Ci which 
contain integrals of the form J dkk l ~ 2 .... The detailed 

expressions for the Cf } ' s are given in Appendix [Dj Here 

we just note that the term C\ represents the lensing 
contribution. As we shall see, it dominates for sufficiently 
high rcdshift and sufficiently large I. Another important 
contribution is C\ which contains the peculiar velocity 
of the emitter, the Doppler term. (It also includes other 
contributions which are, however, always subdominant.) 



The results of this section allow the determination of 
the luminosity distance for a given initial spectrum P-q, (k) 
and given transfer function T k (r)). The transfer function, 
the conformal time rj{z), as well as the conformal Hub- 
ble parameter H(z) depend crucially on the cosmological 
parameters. In a forthcoming paper |l3| we will present 
a code to determine the luminosity distance power spec- 
trum numerically and discuss its dependence on cosmo- 
logical parameters. In this work, were we mainly want 
to present the method, we approximatively calculate the 
power spectrum for a simple case to gain some intuition 
about the order of magnitude of the different terms. 



V. RESULTS FOR A PURE CDM UNIVERSE 

In this section we approximate the luminosity distance 
power spectrum semi-analytically for the simple case of a 
cold dark matter (CDM) universe without cosmological 
constant, O m = 1, ft\ = 0. We assume a scale invariant 
spectrum of initial fluctuations, 



Py(k) = A(kr) y 1 ^ 1 = A, 



n = 1 . 



(76) 



The amplitude A is known from the Wilkinson Microwave 
Anisotropy Probe (WMAP) experiment, A ~ 10~ 10 0. 

In the radiation dominated past of the universe, the 
Bardeen potential is constant on super horizon scales, 
h] < 1 and oscillates and decays like 1 /a 2 cx 1 /rf on sub- 
horizon scales. During matter domination, the Bardeen 
potential is constant [H| ■ To take this gross behavior into 
account, we approximate the transfer function during the 
matter era by 



T k (v)T k (ri') 



Tl 



1 



1 + p(k Vcq y 



(77) 



where ?7 C q denotes the value of conformal time at matter 
and radiation equality. Comparing this rather crude ap- 
proximation with the numerical one, which can be found 
e.g. in Dodelson's book [lij|. we find /3 ~ 3 x 10~ 4 . In 
addition, there is a log-correction which comes from the 
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logarithmic growth of matter perturbations during the 
radiation era. We shall take it into account only for the 

(5) 

dominant term C\ . Furthermore, we use that during 
the matter dominated era 47rGa 2 (p + p) = |(a/a) 2 = 
i(2/77) 2 = 6/r, 2 . 

To determine the power spectrum, we have to perform 
integrals over time of the form 

rv 

1(f) = / dr)'f(r)')je(k(r) - if)) 

= \J x S dx'f( m -x'/k)j t (x') , (78) 

where we have introduced x = k(rio — v)- The spherical 
Bessel function of order £ is peaked at x ~ £. For values 
much smaller than £ it is suppressed like (x/£) and for 
values much larger that £ it oscillates and decays like 1/x. 
In our crude approximation, we neglect contributions to 
this integral from outside the first peak and approximate 
the integral over the first peak by the value of / at x = I 
multiplied by the area under the peak. This gives 

/(/)^Wr?o-|Wfc-— l — W— "*) - 

« V k J \ m - vs ) \m - r) ) 

(79) 

where le is the area under the first peak of the Bessel 
function ji and 9 denotes the Heaviside function, 9(x) = 
0, if x < and 0(x) = 1, if x > 0. Numerically we have 
found ij ~ 1.58/£. Most of the resulting integrals over 
k can either be obtained analytically in terms of hyper- 
geometric functions |15| | or they can be approximated by 
the same method. Finally, one fc-integral contributing to 
the Doppler term C\ has to be performed numerically. 
More details are given in Appendix lEl 

We have tested our approximations by comparing them 
with the numerical result and have found that we nearly 
always overestimate the numerical result, but never by 
more than a factor of 2. The approximations are quite 
bad at low i < 5, but become reasonable later. A fully 
numerical evaluation as we shall perform it in |13| . will 
probably give a somewhat smaller result but not by more 
that a factor of 2 to 4. Here, we are not so much inter- 
ested in numerical accuracy as in qualitative features of 
the different contributions to the power spectrum. 

In Figs. E| to El we show £(£ + l)cf\z, z) for different 

values of z. For i ~ 10, the lensing contribution uf' 
always dominates if z > 0.2. It is interesting to note that 
the different contributions do not scale in the same way 
with £. Only and are scale-invariant with 

£(£ + l)cf ] ~ 10- 10 (80) 
£{£ + l)CP ~ -10- 10 . (81) 

The other contributions grow up to a redshift dependent 
maximum (minimum) from where they decay. They may 
become scale invariant at higher £, but until I = 300 
the scale invariant piece is only clearly visible for z = 
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FIG. 5: The contribution C e . We choose the same colors 
(linestyles) like in Fig. 21 The contributions for z = 0.1,0.5 
and 1 are negative while those for z — 2, 4 pass through 
at low £, visible as a spike. This may well be due to our 
approximative treatment. 



0.1. Higher values of z have their maximum contribution 
at higher £ and have not decayed into a scale invariant 
behavior until £ — 300. The lensing contribution C*- 5 ' 
even just grows. For z = 0.1 it does seem to reach a scale 
invariant plateau, for z = 0.5 it seems just to reach the 
turns over around £ = 300. For values z > 0.5 shown in 
Fig. El the spectrum is simply growing and has not yet 
reached the turn over until £ = 300. 

The most surprising result is the high amplitude of the 
lensing term . Let us discuss this term in more detail. 
After performing the time integrals as outlined above, an 
integral J dkk _1 T^ from £/(rjo — f]s) to infinity is left. If 
we neglect the log in the transfer function, this amounts 
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FIG. 6: The contribution £(£ + l)Cf\z,z)/(2n), without 
the numerical part, for z — 4,2,1,0.5 and 0.1 (from top to 
bottom). Note that here we have chosen linear as opposed to 
a log representation. 
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FIG. 7: The Doppler contribution of 1(1 + l)Cf ) (z, z)/(2tt) 
which has been determined numerically for z = 0.1,0.5,1,2 
and 4 (from top to bottom). Our numerical code is stable 
only for I ~ 80 and we therefore plot only this part of the 
curve. 



to 



dk. 



-Ti 



lo g(#7 

{■no~r)s) 4 



if 



Vo-Vs < /3 1/4 »)c q 



else. 



Together with the factor If x £ 2 from the time integra- 
tions, we obtain a £~ l behavior of 1(1 + l)^ 5 '' at large 
£, which is not seen in Fig. [5J However, when taking 
into account also the log correction, the correct ampli- 
tude and scaling with £ can be estimated in this way (for 
more details see Appendix lD|l . 

This dominant term comes actually from the second 



FIG. 8: The contribution 1(1 + l)C f (4) (z, z)/(2tt) for z 
4, 2, 1, 0.5 and 0.1 (from top to bottom). 




10" 10 t=L 



FIG. 9: The lensing contribution i(l + l)Cf } (z, z)/(2ir) for 
z = 4, 2, 1, 0.5 and 0.1 (from top to bottom). For clarity, we 
have again chosen a log representation in this graph. 



derivatives of '5, hence from the Riemann tensor which 
describes the tidal force field, i.e. geodesic deviations. 

If the fc-integral would not be decaying, 1(1 + 1)C^ 5 ' 
would be growing like ~ £? . But the integrand becomes 
small for fluctuations with wave number smaller than 
about fc cq = l/^VSfeq). Therefore 1(1 + l)cf' has a 



(broad) maximum £ 



q 

max 



k cq (va ~Vs)- Hence £ max is 



increasing with the source redshift. For Zs — 1, hence 
Vo — Vs — 0-3vo — 3077 cq we find £ max ~ 250. The general 
expression for a matter dominated universe is 



e max (z s ) 760 x 



(82) 



Our first important finding is that the tidal force field, 

(5) 

represented by C) totally dominates the final result for 
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Doppler part of l(i + l)Cf \z, z) / (2ty) is important, which 

we have computed only for £ ~ 80. For clarity, we have again 
chosen a log representation in this graph. 
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for z — 0.1 are shown. For this low redshift they are all of the 
same order of magnitude. For low £'s our approximations are 
not trustable, they even lead to negative values for C e ° for 
£ < 3. 



redshifts z$ ~ zgi ~ 0.2. In a numerical treatment, 
where we want to reach a 1% level accuracy, it is suffi- 
cient to consider only for redshifts z$ ~ %s> ~ 0.5. 
Secondly, naively one would expect a result of the or- 
der of (* 2 ) ~ A ~ 10 -10 , but we found nearly 10 -5 for 
supernovae with redshift zs ~ 2. This comes from the 
fact that in the time integral for , the fluctuation is 
multiplied by the conformal distance T) — Tjs- A small 



angular deviation at 77 builds up to a large deviation at 
r]s if the distance is large. Furthermore, we deal with an 
integrated effect where even if the deviation from each 
fluctuation is similar, more small fluctuations pile up on 
the way from the supernovae into the telescope. Even if 
these are uncorrelated, we still gain a factor y/N by pil- 
ing them up. These arguments are somewhat simplistic, 
but they explain, why the term with most time integrals 
and with the factor (77 — 775) dominates. 
In Fig. ^5] we show the sum 



£(£+!) 



1 

2^ 



For z > 0.1, the total results are indistinguishable from 

(5) 

C\ alone. Only for z — 0.1 all terms contribute, espe- 

(3) 

cially the numerical part of C\ dominates. We plot this 
line only until £ = 80 since we have no reliable results 
on the numerical contribution to C\ for higher values 
of I. The different contributions to Ci for z = z' = 0.1 
are shown in more detail in Fig. 1111 

It is also interesting to study the behavior of Ct(z,z') 
for fixed z' as a function of z, and for fixed z ^ z' as a 
function of £. We show this behavior in Figs. 1121 to 1141 
Somewhat surprisingly Ci(z, z') shows no peak at z = z' . 
It is therefore not problematic to include relatively large 
bins Az in a study of Ci(z, z). 
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e yz, z')/(2ir) is shown as 
0.5 and I = 200, 100, 50, 10 



FIG. 12 

a function of z with z fixed to z 
and 2 (from top to bottom). Above £ ~ 50, the ^-dependence 
of the result becomes weak as expected. For z > 0.1 this 
represents actually also the total contribution to Ce. 



VI. CONCLUSIONS AND OUTLOOK 

In this work we have determined the correlation func- 
tion of the luminosity distance fluctuations. We have 
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Again, for z 7^ 0.1 this result is equivalent to the full Ce ■ 



larger than the CMB anisotropy power spectrum! But 
nevertheless, the fluctuations obtained within linear the- 
ory are still much smaller than 1. We have also seen that 
small scale fluctuations do not significantly contribute to 
the Ce's for low £'s. i.e. on large scales. This indicates 
that they cannot change the observed cIl(z) by factors 
of order unity, which would be needed to mimic acceler- 
ated expansion in a matter dominated universe. Also the 
variance, i.e., the typical deviation of a given luminosity 
distance di(n, z) from the mean, which is dominated by 
small scale fluctuations (the lensing contribution) is 

d L {z)- 2 {d L {n, z) 2 ) = -L Y (21 + l)C t ~ 1(T 5 « 1 . 



Our findings thus indicate that the explanation of ac- 
celerated expansion put forward in 6] is probably not 
realized. Of course we have not taken into account the 
change of the transfer function due to nonlinearitics. To 
determine this effect more precisely we would have to 
take into account the non-linearities, especially in the in- 
tegral for cf^ . 

We suggest that the newly derived luminosity distance 
power spectrum given by the Ci{zg, z<?') can be used as a 
new observational tool to determine cosmological param- 
eters. For 1% accuracy of the fluctuations at zs ~ 0.5, 
only C\ has to be taken into account and therefore the 
numerical complexity of the problem seems to be quite 
moderate. In a future paper [13J we shall investigate the 
possibilities to measure Ci(zs, Zs>) with the supernovae 
searches which are presently under way or in planning. 
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found that at redshifts z > 0.2, the result is dominated 
entirely by the 'lensing term' (| A\I/| 2 ) which is propor- 
tional to the density fluctuation. Geometrically it comes 
from the term A? = B?^ vi ^k v i.e. the Riemann tensor. 
Hence this contribution is due to the tidal force field. 
We have seen that it is dominated by fluctuations of the 
size A ~ 77 C q which enter the horizon at matter radiation 
equality. These fluctuations have not been damped dur- 
ing the radiation era, but they are the smallest and there- 
fore the most numerous which have not suffered damping. 
Their effect can therefore add up most along the path of 
the photon. 

We have found that within linear perturbation theory, 
the di-power spectrum is nearly 5 orders of magnitude 



APPENDIX A: CHRISTOFFEL SYMBOLS AND 
THE RIEMANN TENSOR OF SCALAR 
PERTURBATIONS IN NON-EXPANDING 
SPACETIME 

Here we write down the Christoffcl symbols and Rie- 
mann tensor for the metric 

g^dx^dx" = -(1 + 2*)d?y 2 + (1 - 2^)j lj dx l dx 3 

to first order in the gravitational potential 5*. 

r° = * (Ai) 

It = (A2) 
r*oo = (A3) 
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T J oi = -<5f* (A4) R l m = -V'Vj*-^ (All) 



= -7tf* (A5) ff 0jm = iJV,** - J^Vj* (A12) 

r} m - -<5*a m *-^d J * + 7 Jm c> l <I' (A6) ff j0m = - < 5; n V J * + 7jm V J * (A13) 

R Z jmn = - KVjVm - #mVj V„+ 

7 J mV i V„ - 7 J „V i V m ] * . (A14) 

R°ooj = (A7) 
R° 0ij = (A8) 

R iOj — VjVj^ 7ij v I' (A9) Here denotes the covariant derivative w.r.t the metric 

i?°ym = 7ijV TO * - 7™V 3 t (A10) 7y. 

APPENDIX B: THE DERIVATION OF EQ. (J59J) 

We first re-introduce the velocity of the source vg and we collect all terms which contain spatial derivatives of the 
form n l di^ at the end. This brings l(55|l into the form (we dismiss the tilde in this appendix) 

d L (z Sl n) = (1 + zs)|(??o - Vs) + (*o - vo • n) - (??o - Jjsi^ 1 ) *s - (?7o - ?7s - W^ 1 ) v s • n 

rvo rVo rv 

+2 / dr/*- / d?7 / dr)'(r)' - r) S )V 2 ^ 

/■i)o rvo rv rVo rv ~) 

/ d-qn ■ W + 2 / / d?7'n • V* + / dr? / dr?'(r?' - nsWn? didj^ \ . (Bl) 

J?7 S Jr/s Jr)S J*]s Jrjs ! 



2 



Now we use 



— = * + n • V* 

OX] 

to convert all derivatives of the form n • into time derivatives. This leads to 



d L {z Sl n) = (1 + z s )| (770 - T) S ) - ^ (#0 + v • n) + (-2(770 - Vs) + Ug 1 ) * s - (vo - Vs - Ug 1 ) v s • n 

rvo rvo rv 

+(r)o ~ Vs)^o +2 dr/^ - dr) \ , 

Jr/s Jvs Jvs 



dr) (n - Vs)^ ^ 



Us 



VO fVO fVO fV 



I drj%-2 I ^77(77- 775)*+ / dr) / dr)' (r)' - r) S )V } . (B2) 

'Vs •'vs Jvs Jvs 



Via integration by parts we can now convert the double integrals over time into single integrals. For this we use that 

rvo 
'VS 



for a regular function 7(77) integrating by parts § rto dr)(r) — r)s) 2 f(rj) gives 



vo rv rvo 
dr) I dr)'(r)' -r) S )f(r)') = / dr)(r) - r) S )(r)o - v)f(v) ■ 

Vs Jvs Jvs 

Using this in the two double integrals above we obtain Eq. I|59|l . 

APPENDIX C: THE POWER SPECTRUM ^/ x s = 1 f d 3 ke ikx ^ Ck) (C2) 

(2tt) 3 J 

in ,r • , r , • The time evolution of the Bardeen potential is given by 

We use the Courier transform convention \ . ° J 

the transfer function, w(k, 77) = Tk(r))^/(K), which is nor- 
malized such that ^(k, 770) — > \&(k) for k — > 0. Since 
_ ik . x //-is the Bardeen potential is constant on very large scales, 

this identifies ^(k) also with the Bardeen potential right 



*(k) = J d 3 a;e- lk x *(x) , (CI) 
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after inflation. The correlation function 

C*(|x-y|) = <*(x)*(y)) 
depends only on the distance |x — y|, so that we obtain 

(*(k,7?)**(k',7/)) = 

T k ( V )T k ,( V ') yd 3 xd 3 2/C*(|x-y|)e- ik x+lk ' y 

= T fc (7 ? )T fe ,(7 ? ')^ 3 J P*(fc)(2^) 3 <5 3 (k- k') , (C3) 
where we have introduced the power spectrum 



P*(fc) = k 3 J d 3 zC*(z)e^ k - : 



(C4) 



It is easy to verify that this definition is consistent with 
the one given in Eq. |(SSJ. 

Standard inflationary scenarios give ~ A(krj ) n ~ 1 
with n ~ 1. From WMAP and other measurements of 
CMB anisotropies we have A ~ 10~ 10 . We first want 
to determine the correlation of the Bardeen potential at 
positions x = xq — n(j]o — v) an d x ' = x o — n '{vo ~ *?')■ 
With the above we have 

(*(7 7 ,X)*(7 7 ',X')> = 
j3 i„ j3 t,/ 



(27T)* 



e -ik-n(?)o-»?) e +'ik'-n'(?)o-»?') 



Using the identity (see e.g. 0) 

e ik.n( TO -,) = ^- (2f + i)^( fc(7?0 _ 7? )) P ^(k . n ) (C6) 



and Eq. (|C3|) we obtain 

(*(r ? ,x)*( ?7 ',x')) = 
l£(2i + l)(2^ + l)^' / fninnW) 



(2 



je(k( V o-v))k'(HVo-v')) J c^VMk • n)PHk • n') = 
-L^(2£ + l)P,(n-n') / ^P^(k)j e (k( Vo - v)) * 

je(Hvo-v')) = 
-C?\z,z')Pt(n.n') , (C7) 



(C5) 



4tt 



where we have used Eq. 11G6II for the last equal sign. Here 
k is the unit vector in direction k and dQ^ denotes the 
integral over the sphere of k-directions. 

In the same way one derives Eqs. (|67|l to 1)74(1 . Each 
factor in • k can be written as a derivative w.r.t 770 — ry of 
the exponential and therefore replaces je{k(vo — v)) by 
—kj' e (k(r]o ~v))- The Laplacian simply corresponds to a 
factor — k 2 . 



APPENDIX D: DETAILS FOR THE POWER SPECTRUM 

In this appendix we write down in detail the expressions for the Cjp y s used in this paper. 

As mentioned in Section IIVI the power spectrum of the luminosity distance can be split in five different parts 
containing ^-integrals of different powers, 

contains the integrals of the form j represents the redshift and parts of the integrated contributions. 
C^P contains the integrals of the form J dk, represents the correlation of the Doppler term with the terms in . 

C'P contains the integrals of the form J dk ■ k, represents the Doppler term and some (subdominant) integrated 
terms. 

Ci 4 ^ contains the integrals of the form J dk ■ k 2 , is dominated by the correlation of the Doppler term with the 
lensing contribution. 

contains the integrals of the form J dk ■ k 3 , represents the lensing term. 

From Eqs. (|5l)|) to (|74*|) and the expression lt58|) for the luminosity distance we obtain the following expressions for 
the C £ (i) 's 



a 



(i) _ 



2 f v ° ( 1 

/ dr)T k (r))jz(k(r) - 77)) - 1 + — — 

1o-i)si, s V ^s{Vo~Vs) 



Tk(vs)je(Hvo - Vs)) 



Vo - vs> J Vs 



dr)Tk(r))jt(k(vo - v)) 



1 



1 



Ti-s'ivo ~ Vs>) 



Tk(r)s')je(k(rio ~Vs>)) 



(Dl) 



C, 



(2) 



dkP^,(k) 



3H S 



Ksivo - Vs) 



(T k (vs) + H^Tfcfos)) j' e (k(Vo ~ Vs)) 
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1 f v ° 1 
■ — — ; / dr]Tk(r])j' e (k(r)Q — rj)) H 

^siVO - Vs) Jr,s % - VS Jris 



Vo - Vs' J Va 



dr}T k {ri)ji{k('no -??))- 1 + 



77s <^> ?7s' 



vo rv 

dv / dv'nWftikivo - V')) 
J vs 



Tk{vs')je(Hvo -VS')) 



(D2) 



Cf> = ~ I dk fcP*(fc) 

7T 



1 



1 



3ft s V Ksivo - Vs) 



(T k (r,s) + nffkins)) j'Mvo - vs)) 



i p i 

■—— / dr)Tk(v)fi(Hvo ~ V)) H , 



vo rv 

d V / dv'T k (v')j'MVQ-v')) 
J vs 



m S ' 



i 



fts'ivo ~ Vs>) 

1 ffo 



(TfcfasO +Hs}T k {n s .)) j'AMvo - VS')) 



+ 



Ks'ivo ~ Vs') Jn s , 
2 r^o 



dvTk(v)fAHvo - v)) + 



1 /■»7o rv 

/ dv dv'T k (v')j'AHvo-v')) 

Vo ~ Vs' J Vs , J Vs , 



n{vo - Vs) 

2 rvo 



Vs J VS 



dk fcP*(fc) I dr\ J dv'{v' ~ Vs)T k {v'){jt{k{vo - v')) +fe(HVo ~ v')) 

1 



Vo - VS' J Vs , 

+ Vs VS' 



dvT k (v)je(k(va -'?))- ! + 



T-Ls'ivo ~ Vs') 



T k {vs')ji{Kvo -VS')) 



(D3) 



C, 



(4) 



7T r?o - ?7S' 



& k 2 P^{k) 



3H S V ^s(?7o - ?7.s} 



(T fc (%) + n^t^s)) j'Mvo - vs)) 



■— -/ d V T k ( V )j' e (k(vo - v)) + / dv d V 'T k (v')j'e(k(vo-v')) 

Ks{vo ~ Vs) J vs Vo ~ Vs J^s Jvs 



Vo rV 



Vs' J Vs 



dV / dv'W ~ VS')T k (v')(ji(k(vo - v'))+j"(Hvo ~ v'))) + 



VS ^ VS' 



(D4) 



C, 



(5) 



dk k 3 P^(k) 



7T {vo - Vs)(vo ~ VS') 

dv j dv{v ~ Vs)T k {v)\3i{k{vo - v))+je(Hvo - v'))) 



Vs J vs 
Vo rv 



Vs' J Vs 



dr, / dv'iv' ~ Vs>)T k (v')(h(k(vo - v'))+3"{Hvo - v')) 



(D5) 



APPENDIX E: INTEGRALS AND APPROXIMATIONS 

Here we make full use of the relatively crude approximation l|79f) 

dxf(x)j e (x) ~ hf{i)6{x 2 - £)9(£ - x x ) , 



(El) 



where 9 denotes the Heaviside function, 6{x) = 1 if x > and 8(x) = else. Hence we neglect contributions to the 
integral which do not come from the region of the first peak of the Bessel function. This procedure is very useful to 
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estimate the result, but cannot be trusted better than within a factor of about 2. We have tested it with numerical 
examples |l6j . A more detailed numerical treatment will be presented elsewhere [lot . Furthermore, we assume a 
scale-invariant spectrum with = A ~ 1CP 10 . We also use the fact that in a matter dominated universe the transfer 
function does not depend on time and can be taken outside the time-integrals. 

We define bs = ^ = x y 1 \_ zs , x s = k(r) — rj S ) and as = (3 f ) ■ Note that xs' = ^ xs- In terms of these 
variables, the transfer function becomes 

T 2 (x s ) = — - 4 , (E2) 
1 + asx s 

(5) 

except for the C\ , where we have to take into account the log-correction. 

1. C« 



r {i)< x 2A (2-6 s )(2-6 s Q f°°dx Sr ^ f , 

c * {zs > zs,) = 4(1 -h s) y — r (^)^(^)^fe) 

+ 4(1 - 6 S )^°° ^r-T 2 (x s ) Qf " dxj t (x)j Qf *' dxMx)^ 

- (2- b s >) £ ^r T2 (*s) (jf* dxjtixfj ■ u{x s .) + b s ^bs^. (E3) 

For the first term, the integral converges without the transfer function, we may therefore neglect it and perform the 
integral analytically. For the second and the third term, we use the approximation IjEljl. Assuming that zs < zs> (if 
not, we reverse zs and zs> in the formula), we obtain 



C «fc, Zq A *±[^X-bs_ rdx s 1 I± 2-bs 1 

1 {ZS,ZS ' ~ tt y le 1 - b s > h 4 1 + «sx% Pl-b s ,l + l^a s 



-(2 - b s )(2 - b s >) y °\ " ,J F [e,£+l;2£ + 2 



16 T(l + 3/2) v " "°' J (2 - b s - b s >) 2e V (2-bs-bs') 2 

(E4) 

Here F denotes the hyper-geometric function and T is the T-function. We use the notation and normalization of . 



2. G, 



r (2), s -2A J (2-3fe s )(2-b y ) rdx s 2 

C) '(Z S ,Z S >) = —r 777- 7-7 / T {xs)ji(xs)ji(xs) 

vr(l - OsO { 2(1 - 6 S ) 7 x s 



12(1 -&s) 2 

- 4(1 - b s ) J™ ^rT 2 (x s ) [J*' dx Je (x^j f^J* 3 ' dx Je (x] 

+ (6 - 7b s >) I" ^T 2 (x s ) ( H dxMx)) ■ j t (x s >) 
Jo x s \Jo J 
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Here again, the terms which contain only an integral over xs can be calculated analytically when we neglect the 
decay of the transfer function. For the other terms we use the approximation (|E1|) . 



Cf\z s ,z s ,) = 



-2A 



7r(l-6 s )(l-6s0 

G i 



- 8/|(l - b s ) 



dx s 



£ 2 1 + a s £ 4 

h b s (2-Zb s ) 



(6-7b s )(l-b s ) + 
1 



g x| 1+ asXg 
Ab s (2 - 36 5 ) a s l A 



£ 2 (2I + 1) 



1 + a s £ 4 



-(£ + iyi e+1 e [e-(e + i) 



l-b s 



1 + a s £ 4 



•-{1-1) 



l-b s 
l-b s > 



6 S /(2-36sO 1 



£ 2 (2£+l) 3 l + a s >£ 4 

-{£+\fh +1 9 (£-(£+!) 



l-b s 
£{£-l)h-i 
l-bs 



£-{£-!) 



bs> 



l-bs 



l-bs 

(4(1 -b s - b s >) + 3b s b s ,)— r(£ + 3/2) {2 _ bs _ bs , r 
bs(2-3b s ){2-b S ')\V^ r(£+l) (l-bs) 1 - 1 



F 



12(2£ + 1) 



2 T(£+l/2) (l-bs'Y 



Fli, -1/2; £+1/2; 



l;2£+2; 

(l-bs? 
(1-bs'f 



4(l-b s )(l-b s ,) 
(2-b s -bs-) 2 



na r(£ + 2) (i-&s)* +1 



F [£+ l,l/2;£+5/2: 



(1 



4 r(£ + 5/2)(l-& S ')^ 
6 S ,(2 - 36 S 0(2 - b s ) r(£ + 1) (1 - 



(l-6sO s 



12(2£+1) 



4 r(£ + 3/2) (l-b S 'Y +1 ' 



F £,l/2;£ + 3/2: 



(1 



(1-bs'f 



A r(£ + 2) (i-6 s )* 

2 r(£ + 3/2)(l-6 s 0' +1 ' 



F £ + 1,-1/2; £ + 3/2; 



(1 



6s) 2 



(1 - 6 £ 



(E6) 



3. C. 



(3) 



Cj 3) (zs,2s0 



&s(2-3b s )(2-3bsQ 

6(1-65) JO 
6565/(2 -36s)(2-36 s Q 
36(1 -6s) 2 

4(1 - 6 S ) 2 r ^fT 2 (x s ) ( dxj £ (x) 
Jo x s \Jo 

36 S '(1 - 6 S ) r ^T 2 (x s ) ( r dxj e (x)) ■ jt(x s >) + b s «*• 6 £ 

JO ^5 VJO / 



dxsT 2 (x s )j'i(xs)ji(xs') + 6 S 44> 6 S ' 
dx s x s T 2 (xs)j'i(xs)j'i(xs>) 



dxj e (x) 



bs<(2-3b s >) f x dx s 
3 Jo 



T 2 (x s ) / dxjt(x) ■ jKxs') + 6 s ^ 6 s 



+ 2(1 - b s ) 2 J™ ^-T 2 (x s ) Qf" dx£ dx'(x s - x')j e (x'f) ■ (^£ S ' dxj t (x)j + b s & b s > 
(2-b s ,)(l-b s ) J°° ^T 2 (x s ) (J** dx£ dx'(x s - x'Mx'fj ■ jtips') + bs «*• b 8 \ (E7) 
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Here, it is not possible to neglect the transfer function in the third integral, because for z$ — zs> the integral does 
not converge without T 2 (xs). We therefore have to calculate the third (the Doppler term) term numerically. 



9 A I f c 

C ^ ZS ^ = 271(1 - bs) J< 

i} i 



dx s 



+ 



+ 



+ 



P 1 + a s £ 4 

It b s (2-3b s ) 



3bs(l-bs)-£ 2 {l 
1 



b i)(bs-bs^ + bs(bs- 2/3)^^ 



£ 2 {2£+l) 



1 + a s £ 4 



-l)I e -iO[l-(£-l) 



l-b s 
l-b s > 



-{£ + l) 2 h +1 e[£ -{£ + !) 



l-b s > 



h b s >(2-3b s >) 



£ 2 (2£+l) 3 l + a s ,£ 4 

-{£+lf h +l 9 [£-{£+!)'- _ h 



£{£-l)h-i6 [£-{£- I) 



1 - bs> 
l-b s 



1 - be, 



+ (2 - -(6 S + b s >) + 8b s b S ')- 



T(£) (l-bsY(l-b S 'Y 



+ 



2 V ~" 1 ' ~~"~" / 4 r(£ + 3/2) (2-6 s -&5<r 
6 s (2-3fe s )(2-36 s 0rV^ r(*+l) (1 - bsY' 1 



F e,e+l;2£ + 2; 



4(l-6 g )(l-by) 
(2-6 s -6s0 2 



6(2^+1) 



2 r(^+l/2) (l-6s0 £ 



F l,-l/2;l+\/2\ 



* r(, + 2) < 1 - ta '' t \p^ + i,W + «A (1 -'* ) ' 



+ 



4 r(^ + 5/2)(l-6 s 0' +2 V >'"""' '(1-6 S ') 2 
b S '(2-3b s >)(2-3b s ) [V^ r(£+l) (1-6 S ) 



6(2^+1) 



4 r(£ + 3/2) (i_6 s ,)m- 



F £,1/2; £ + 3/2; 



(l-bs') 2 



(1-M 5 



' 2 r(£ + 3/2)(l-6 s ,)' +1 V 7 7 '(l-^') 2 



+ 



6 s 6 s ,(2-36 s )(2-36s') 



dx s x s T 2 (x s )j' e {x s )j' e {x s > ) 



36(1 - 6 S ) 2 

The last term in this sum is determined by numerical integration over x s . 



(E8) 



a 



(4) 



c?w) - , (1 - M(i - M { ( 2 - 3 ^)( i - »»> r ^(^(X5^ y ) 

+ ' s(2 6 ( ?- * g ° jf ^(^)i/(^)^(^0 

+ 4(1 - 6 S ) 2 ^°° ^rT 2 (x s ) (jf " dxj t {x)j f^fj' dxj t {x)J 
- 2(3-4bs>)(l-b s )J~ ^-T 2 {x s )(^£ dxuix^j .j t (x s >) 

- 2(1 - bs) 2 £ ^T 2 ( x s) (jf * dxje(xf) ■ dx £' dx'(x s , - x')j e (x') 

+ (2 - 3b s )(l - 6s) I" ^T^xs) ( r' dx H' dx'(x s , - x')j e (x')) ■ j t (x s ) 

JO X S \J0 Jx / 
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+ 



bs(2-3b s ) 



r ^t\x S ) ( r dx r dx'{ XSI x 1 ) 

Jo X S \Jo Jx 



)j t {x')\ ■ j[{xs) 



+ b s b s > . 



(E9) 



C { e 4 \z s ,z sl ) 



-2A 



tt(1-6s)(1-&s') 



t 2 i 
1 i ■ 



' 2/?(1 - 45)2 f ^ TT^J ( 2 < 2 + 1 '> + ^Ts ib ° + »* - 2) 



Aa s £ 4 



l + as£4 [ ~ 2(1 - 6s)(3 - 4fe s ) + £ 2 (2 - 3fo s )(6 s - b s ,) + b s (b s - 2/3) i + ^ 4 



+ 



+ 



// 65(2-365) 1 



£ 2 (2£+l) 3 l + a s £ 4 



-{£+l) 2 h +1 e[£ -{£+!) 



l)I e -i9 [£-{£- I) 



1-bf 



l-b s > 



±-b s > 



h 6 S '(2-36 S /) 1 



£ 2 (2€ + 1) 3 l + a s >£ A 



-(£+l) 2 h +1 9[£ -{£+!) 



'-I)h-i0 [£-{£- I) 



l-b s > 



l-b s 

h£ b s (l-b s ) 3 (2-'Sb s ) l £I l _ 1 b s >-l + £(bs-b s >) 



2£+l 6 

l-bs>+t(bs-b s > 



it+i 



£-l(l-b s )* + (£-l)%i q 
l-b s 



l-b s > 
l-bs 



t-i-e 



l-b s > 



(i-b S r + (£+im q 



£+1 



l-b s > 



+ b s & b s > 



+ 



4 r(£ + 3/2) {2-b s -b s >) 2e 
6s(2-36s)(l-Mrv^ r(* + l) (1-&S)'" 1 



(2-fe s -6 s 2 



6(2^+1) 
/5F r(£ + 2) (l-6 s )' +1 



2 r ( / + i/2) ci -?,,,)' 1/V + 1/2; (i-&s/) 2 



F 1,1/2; £+5/2: 



+ 



4 r(£ + 5/2)(l-& s 0" +:r 
6s'(2-36fir0(l-6s)rv^ r(*+l) (1 - 6 S ) 



(i-M 2 

(1-M 5 



6(2^+1) 



4 r(£ + 3/2) (i-b s ,)W 



F £, 1/2; / + 3/2; 



(i-M 2 
(i-M 2 



^ r ^ + 2 ^ ^ " F f/ + 1,-1/2;/ + 3/2; 



2 r(£ + 3/2) (l-6s0" +1 " 



(1-M 2 



(E10) 



5. C. 



(s) 



C| 5) (zs^s') 



2A 



tt(1 - b s )(l - b s >) 



(l-6s)(l-6s0 



/'°° ctes 
Jo 
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dxj e (x) 



- 2(1 - b s ) 2 j°° ^T 2 (x s ) ( r' dxje(x)) ■ jt(x s ) + b s ^ b s > 
Jo x s \Jo / 

- 2(1 - b s ) 2 ^-T 2 (x s ) (jf *' dxj e (x)j QT* dx£ S dx'(x s - x')je(x')^ +b s ^ b s > 
+ (l-bs)(l-b s -) r ^T 2 (x s ) ( r dx n dx'(xs-x')j e (x'))je(xs') + bs^bs< 

Jo X S \Jo Jx / 
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+ (1 - b s ) 2 J™ ^T-T 2 (x s ) Qf * dx£ dx'{x s - x')j e (x'?) (jf *' dx£' dx'{x s , - x')j t (x')j J . 

(Ell) 

The first term is dominated on large scale and we may thus set T = 1 so that it can be integrated analytically. For the 
other terms we use again the approximation (|iill|) for the integrals dx or dx' . The biggest contribution then comes from 
the last term where we have to perform two double integrals dxdx' , which result in lf(2 + £ 2 — £xs)(2 + £ 2 — Ixs') oc £ 3 . 
In this term, which becomes large for large £ or lar ge x s, we take into account the log-correction to the transfer function 
for better accuracy. From the expression in Ref. 14] and our definitions we find: 



T 2 {x s ) = . (E12) 

1 I usi s 

in 2 (i+ -?i 1 "; 4 ^T 



Using our approximation IIE1I) . we obtain 



<#W) = ^^U 2 (l-M 2 f f 1 —< (2 + £ 2 -£ Xs){ 2 + £ 2 -£x s ,) 



(2- {2 + ( 2 )b s + £ 2 b s ) 



If l-&s 

£ 2 1 + a s £ A 

^ r(£) (l-^l-M^^ / + 1 .^ + 2 .4(l-y( 1 -6 f )- _ 



where 



4 r(^ + 3/2) (2-6 s -M 2£ V (2-65-&S') 2 



«S=/3(t^) 7 1 T- (E-L4) 

■■■ 1 - S/, , 7.8-10- A " X 



In* 1 + 



1-fcs 



(5) 

The remaining integral represents by far the largest contribution to C] . For sources with equal rcdshifts z$ — 
Zs> = z, the spectrum ci 5 \z, z) grows until a s £ A ~ 1 and decays for larger £. Neglecting the log correction we have 

(.4 
/3 1 / 4 j = ^max- Hence C^ 5 '' grows roughly until £ max and decays afterwards. With & oq = (?7oq/?7o) — 0.01 

we obtain 



I ~ 760 VZS + 1 ~ 1 

For a crude order of magnitude estimate, we first neglect the log corrections. For £ <C £ max the integral is dominated 
by the region xs < W and we may simply integrate until xs — W, neglecting the x A s decay of the transfer function. 
In the opposite region, if £ 3> £ max , we may neglect the 1 in the denominator of the integral. An interpolation between 
this two asymptotic regimes gives 

g(5) f 2^ [ln(V) + i if £<l max 

Cl {ZS ' ZS) - - if (E15) 

Since I 2 cx l/£we see that £(f+l)C'l 5 ' ) grows like ^ 3 for small £'s and it decays like l/£ for large £'s. The broad maximum 
is reached roughly at £ max ~ 759 Vi+|£ ^i _ ygg^ _ fr s ) and is of the order of (^4/7r)^f nax . This approximation is, 
however, surprisingly bad. We therefore take into account the log in the transfer function by simple replacing as 
by a$, where £ in the expression for as denotes the lower boundary of the integral. The expression for £ max then 
becomes ^-dependent, 



VMl + 7.8xlO- 4 £/(l-M (1 _ 5sr 



cq 
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For £ < 1.3 x 10 3 (1 — bg) = £g the log can be expanded and £ max /£ behaves like £ 1 / 2 leading to a linear growth of 



£{£ + 1)C^ . Only above £3 it levels off. For zg = 2, the asymptotic regime, where £(£ + l)Cf decays like 1/1 is 
actually only reached at £ ~ 2000, where our approximations (and linear perturbation theory) no longer hold. 

In Fig. 1151 we plot the approximation given in Eq. i|E15(l with £ max given in (|E16|) for zg = zg> = 2 and hence 
£g ~ 540. Actually, to have a better fit with the numerical integral we choose a slightly modified value, namely 

£max — 0.75£max- 
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FIG. 15: The approximation for £(£ + l)Cf> (z, z')/(2n) given in Eq. 1E15II (red, dashed line) is compared with our numerical 
result (black, solid line) for z = 2. 
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